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The appearance of a geometric flow in the conservation law of particle number in classical par¬ 
ticle diffusion and in the conservation law of probability in quantum mechanics is discussed in the 
geometrical environment of a two-dimensional curved surface with thickness e embedded in Rs. In 
such a system with a small thickness e, the usual two-dimensional conservation law does not hold 
and we find an anomaly. The anomalous term is obtained by the expansion of e. We find that this 
term has a Gaussian and mean curvature dependence and can be written as the total divergence 
of some geometric flow. We then have a new conservation law by adding the geometric flow to the 
original one. This fact holds in both classical and quantum mechanics when we confine particles to 
a curved surface with a small thickness. 

PACS numbers: 87.10.-e, 02.40.Hw, 02.40.Ma, 82.40.Ck 


1. Introduction 

The motion of particles on a given curved surface 
is an interesting problem in a wide range fields in 
physics Q-i- The classical diffusion equation and the 
Schrodinger equation on such a manifold is expressed by 
changing the Laplacian to the Laplace-Beltrami operator 
in the equation; however, when the surface has a thick¬ 
ness e, i.e., the configuration space is x , the situ¬ 
ation is not simple [3|- 

Such a 2-1-1 dimensional system can be found as the 
motion of a protein in a lipid-bilayer (cell membrane) 
in classical mechanics [I|. The direct application to the 
quantum mechanical particle is not yet done, however, 
the motion of an electron in the quantum hall device 
may be a good example to test the geometric effect by 
bending the surface of the system. Though it is not for 
the real particle, such a geometrical effect is discussed 
by several authors for the Josephson junction device Q- 
The curvature of junction surface gives the effect to the 
kink motion that is the solution of sine-Gordon equation, 
the effective theory for the phase difference between su¬ 
perconducting electrodes. 

In this paper we consider the conservation law in 
such a pseudo two dimensional geometry and show 
the existence of anomalous term. The case of classical 
mechanics is already discussed in the previous paper 0, 
3, so we just give the result and we devote to discuss 
the quantum mechanical case. 

2. Geometrical Tools 

To make the problem concrete, we first introduce a 
two-dimensional curved manifold E in , and we also 
introduce two similar copies of E denoted E' and E and 
place them on both sides of E at a small distance of e/2. 

Our physical space is between these two surfaces E' 
and E. We set the coordinate system as follows. 

X is a Cartesian coordinate in i? 3 . a: is a Cartesian 
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FIG. 1: curved surface with thickness e 


coordinate that specifies only the points on E. g* is a 
curved coordinate on E, where the small Latin indices 
run from 1 to 2. is the coordinate in R^ nor¬ 
mal to E. Furthermore by using the normal unit vector 
on E at the point {q^,q^), we can identify any 
point between the two surfaces E' and E by the following 
thin-layer approximation Q: 


X{q°,q\q^)=x{q\q^) + q°n{q\q^), (1) 

where —e/2 < q^ < e/2. 

Then we obtain the curvilinear coordinate system be¬ 
tween two surfaces (c R 3 ) by using the coordinate q^ = 
(q°, q^,q^) and the metric (Hereafter, Greek indices 
- ■ ■ run from 0 to 2.) 


dX dX 
~ dqt^ ' dq''' 

Each part of is expressed as follows: 


( 2 ) 


^ r^.dx dn dx dn, , dfi dn 
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where 


dx dx 


( 4 ) 


is the metric (first fundamental tensor) on S. Hereafter, 
the indices are lowered or raised by gij and its 

inverse . We also obtain 


Goi — Gio — 0, Goo — 1- (5) 

We can perform the calculation by using new variables. 
We first define the tangential vector to S by 


G = detGy = g + 2gKq° + g{K'^ + i?)(g°)^ + 0{{q°)^). 

(13) 


The inverse metric of G^ is given as 

Q^J = gV _ + |(2kk®J' - Rg^^){q°f + G(((?°)^). 

(14) 


Furthermore, 

(«;-i)b = |(^^d_^b)^ (15) 

and from this relationship we obtain 


Bk = 



( 6 ) 


Note that n ■ Bk = Q. Then we obtain two relations: 
the Gauss equation 


OB, 

dqi 


—KijU + T^jBk 


and the Weingarten equation 


dn 


= KTBr, 


( 7 ) 

( 8 ) 


^Rg^^ = (16) 

By using the above relations, we can construct the 
diffusion equation and Schrodinger equation in our 
environment. 

3. Effective Schrodinger equation and Geometric Flow 

For the classical diffusion field, the problem has al¬ 
ready been solved and discussed in We just note the 
important result. 


where 

Ty = + 9jgim - dmgij)- 

Kij is a symmetric tensor called the Euler-Schauten 
tensor, or the second fundamental tensor defined by the 
above two equations. Furthermore, the mean curvature 
is given by 

K = g"^Kij, (9) 


Ot oq^ 

f) 1 f) 

= -Vf>(J- + Ji), (17) 

where H, D, and are the diffusion constant, 
~ 2 

D = the two-dimensional effective diffusion field 

and the two-dimensional covariant derivative respec¬ 
tively. The normal diffusion flow J is 


and the Ricci scalar R (Gaussian curvature) is defined 

by 


r = -Dg^^ 


dqi ’ 


(18) 


R/2 = dei{g'‘^Kkj) = det(K®) = ^(^^ ~ k^-k®/). (10) 

Then we have the following formula for the metric of 
curvilinear coordinates in a neighborhood of E: 

Gij = 9ij + 2,q^ Kij KimK^. (11) 

Now we have a total metric tensor such as 

= ( 0 Gy ) ■ 

From the definition of the metric (HB and the Ricci 
scalar (US, we can construct the following geometrical 
quantities: 


and the geometric diffusion flow Jq is 

Jh = -D[ {‘iK^^Ki,-2KK^^f-^-\g^^^ct>^^\ (19) 

The two conditions that the surface is curved, and e 0 
are essential for the existence of geometric flow. 

In quantum mechanics, the basic equation is the 
Schrodinger equation which is written using curvilinear 
coordinates in a three-dimensional space between the two 
surfaces E' and E. 
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where the form of the Laplace-Beltrami operator is 
A(3) = 

and we suppose that V depends on neither t nor (p. 

Starting from this wave function^, we construct the 
effective two-dimensional theory Hi ) 0 ) @ ■ From the nor¬ 
malization condition we obtain 

1 = ^ = y[y ^ 

( 21 ) 

Our effective two-dimensional wave function should 
satisfy 


Fo = (30) 

Fi = /c(i?-y)-i(A(2)«), (31) 

+ ^g"^{diK){djK) + {W,{K"^djK)} 

(32) 

ii = (33) 

A2 = SV.K^'^Kldj. (34) 



Then how can we obtain the dynamical equation for (p 
? To solve this problem, we first define a new variable ip 
as 

^ = (G/5)'/V (24) 


Then our equation for ip is given as follows: 

+g°ii + (g°)2i2]^ + F^, (35) 

where we have omitted 0 {{q^)^) terms for small q^. 

We treat this system by the perturbation method. The 
Hamiltonian can be written as 


with 

\^{q\q^)f = r ^ dq\ (25) 

J-e/2 

Furthermore we suppose that it is possible to separate 
the variables. (Later we prove the possibility of separa¬ 
tion of variable up to second order of e expansion.) 

4’ = V5('7^9^1) x(9°,0: (26) 

^+e/2 

1 = / IxP dq°. (27) 

J-el2 

Then we can construct the equation for tp from the 
Schrodinger equation of tp. To start this program, we 
first construct the Schrodinger equation for ip. This has 
the same form as (1^0)1 except that the Laplace-Beltrami 
operator is changed to the following operator. 

A(3) = (G/5)1/‘‘A(3)(G/5)-iA. (28) 

Using some tools prepared in section 2, this operator 
can be expanded as 


H = Ho+Hi, (36) 

+ dW ^ ^ 

Hi = -^[q°{Vi+A,) + {q°f{V2+A2)]. (38) 


For simplicity we write 

e = 9 °/e, ( 39 ) 

y = {q\q^), ( 40 ) 

h{y,dy) = -| 1 [A( 2 )+Fo] + F, ( 41 ) 

f{y,dy) = -^{V,+A,), ( 42 ) 

Hy^dy) = -^(^ 2 +^ 2 )- ( 43 ) 

( 44 ) 


Furthermore we utilise TA/m = 1 unit in a while. Then 
Hamiltonian can be written as 


A(3) = a( 2) + _|^ + Ho + gVi + (q°) ^2 

+g°ii + (g‘’)'i2 + G((g°)'), (29) 


H = -^dl + h{y,dy) + e^f{y,dy) + e^^^g{y,dy). (45) 
The boundary condition in ^ direction is 


where Fq, Fi, F 2 , Ai, and A 2 are given by 


^(e = ±l/2) = 0. 


(46) 
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We look for a wave function in the form 

^ = '00 + e'01 + e^V’2 H-) 


(47) 


In the same way, we suppose the existence of the eigen 
equation for y direction i.e. 


and energy eigen value 

E = e“^i5(_2) + + i5(o) + ei^(i) + H-• 

(48) 

Putting all expansions into the time independent 
Schrodinger equation, and we obtain the following equa¬ 
tions in each order of e. 

-^'9|V'o = £^(-2 )V'o, (49) 

= ■E^(-2)V'i + ^^(-i)'0o, (50) 

-i|9|V'2 + (iV'o 

= '®(0)V'0 + ^'(-l)'01 + £'(-2)'02, (51) 

= + ^'(o)'01 + £'(-1)V'2 + £'(_2)V’3 j (52) 

-^i9|V'4 + H 2 + ?/Vl + 

= ■E(2)V'0 + E^i)ll)i -I- il(0)V'2 

-f £;(_i)^/> 3 + il(-2)'*/'4- (53) 


h V’niy) = A„ (^„(?/), (59) 

and we suppose the the orthonormal condition for (/?, 

J ^*m{y)‘fn{y)y/g d'^y = s^n- (60) 

From equation (l50l) . we have 

(--i9| - En)4’1 = i?(_i)aoXAr. 

We expand i/'i as 

V'l = (y)xM(C)- 

M 

and by putting this into the previous equation and by 
using dSHl), we obtain 

a^^iyEK — En) = aoE(^_i')6KN- 

when K = N, we have 7?(-i) = 0. And when K ^ N^ 
we obtain = 0. Note that this equation does not 
determine = ai{y). So we have 


The first equation (|4^ is easily solved by using bound¬ 
ary condition (|15)) and we obtain 


XN 


^/ 2 cos{N^T^) A^ = odd, 
•\/2sin(A^7r^) N = even. 


(54) 



with normalization condition 


(55) 


So we have 



\Xn? = 1. 


(56) 


tpo = ao(y)Xiv(0: ^(-2) = Ew. (57) 

Hereafter we utilise the orthonormal condition for y. 


.1/2 


1 - 1/2 


XMiOxNiOd^ = 5 mn- 


(58) 


Note that the boundary condition (H51) is an ideal¬ 
ized condition because it is derived by the infinitely deep 
square well potential. Therefore (l54|l - ([5^ are the ap¬ 
proximated relation by using a deep enough but finite 
confining potential. 


= ai(y)xAr, ^(.i) = 0. (61) 

From equation (ISTI) . (l57l) . (I6T]) we have 

(-ia| - Em)'iI}2 + XNih- E^o))ao = 0. 

We expand '4’2 as 

V'2 = (y)xM(O- 

M 

and by putting this into the previous equation and by 
using (EHl), we obtain 

— En) = —{h — A(o))ao(5ifAr. 

when K = N, we have (h — ill(o))ao = 0- And when 
K ^ N, we obtain = 0. Just like (l6T]) this equation 
does not determine a}q = 02 (y). So we have 

f/'2 = a2{y)xN{0- ao = <7’n, ^^(o) = A„. (62) 

From equation (I52|), dSll), (El]) , El), we have 

(-^9| - EN)'tll3 = Xn{K - E)ai +XN{E(i)ipn - 
We expand i/'s as 

V's = (y)xM(O- 

M 
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By putting this into the previous equation and by using By putting this expansion into the previous equation and 
we obtain by using (1551) . we obtain 


‘‘K 


{Ek — En) = —{h — Xn)ai5KN 


~ ^n) + — A„)a2 


+ SKNE(^i)ipn- < K\^\N > flfn, 


where 


< K\^\N >= J Xk i Xn d^- 
When K ^ N, we obtain 

I^K — 

When K = we have 

-{h - Xn)ai + E(^i)(pn =< -/VICI-/V > fipn = 0. 

Last equality comes from <7V|^|7V>=0. By using the 
expansion 


“1 = 


multiplying y/g^p^, and integrating by y, we obtain from 

dmi) 


0 — bl^\Xn — Xm) + Eti\Snm- 


Then we have 

bl2 = 0 (to n), £ 1 ( 1 ) = 0. 

Note that still remains as undetermined constant, 
however, such term appears as 

■01 = bll'>(pnXN- 

This is the same function as '0o in expansion. So we take 
bn^ = 0 in the following. In total we obtain 

, ^ <mN> ^ 

03 = asXN - 2 ^ E - E 

01 = 0, £1(1) = 0. (63) 

From equation (1551) . (1571) . (ICT|) . (15^ . (1551) . we have 

-^ 0|04 + ha2XN + ^'^WnXN 

= E(^2)‘PnXN + Xna2XN + Ej^lpi- 
We expand 04 as 

04 = (2/)xm(0. 

M 


+ < K\^‘^\N > gipn — dKNE{2)P>n = 0 , 

For K N, we obtain 

For K = N, we obtain just like before 
< 7V|^^|iV > 

02 = - 2^ —;-^- {m\g\n)ipm, 


m^n 




where 


Ei2) =< IVIC'IIV > (n|5l^), 


(TO|g|n) = j Va (t^hgft^nd^y. 


Then we obtain 


04 = a4XAr - 


KyLN 


< K\(^\N> 
Ek — En 


-XK gvn, 


02 = - 2 ^ —;-^- [m\g\n)ipm Xn, 


m^n 




E{ 2 ) = < N\^'^\N > {n\g\n). 


(64) 


In total we obtain the wave function and energy eigen¬ 
value up to O(e^). 


-Ip = {ifn-e 2^ ^- {m\g\n)(pm} 


m^n 

X Xn exjp{—iEt/h), (65) 

E = e~^EN + Xn + e^ < iV|^^|iV > (n|g|n), (66) 

where we brought out the explicit time dependence. 
Note that we have no N number excitation, and separa¬ 
tion of variable method holds up to ^(e^). 

The effective Hamiltonian that leads to (1551) . (1551) is 

Heff = -^0| + hE^< IVie^lIV > g. 

with a solution in the form of 0 = Lp{t) XN(t). Then the 
Schrodinger equation for ip field is 

th^ = [h + e^ <N\e\N>g]ip. 


(67) 
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In the original variable we have 

ih.^ = + Vb + e'^{V2 + A2)}(p + Vip, (68) 

where 


where J and p are three dimensional flow and density re¬ 
spectively. We utilize the cylindrical coordinates (r, 0, z) 
and 

= dr^ + d6^ + dz^, (75) 

= diag(l,r^, 1), ^G = r. (76) 


f^s«»(IV|f>|IV) = d(l-^). (69) 

From this Schrddinger equation, we obtain 

-I- e'^{(p*A2ip — <pA2ip*)} 

= - ifdjp*)} 

= V,(J* + J^), (70) 

where 

J" = (71) 

^ (72) 

Note that the back reaction from ip field to the x field 
occurs in G(e^) (See ((63l) and (l64l) l. 

5. Toy model 

To consider the physical meaning of geometric flow, 
let us show one simple example: bending ribbon with 
thickness e as seen in figure 2. 



- 'Ipdrlp*), 

Zmi 

(77) 

{tp*dgllj Ipdg'tp*), 

Zmr^i 

(78) 


(79) 


Then the conservation law gives 

deJ^ + drr + d,r + -r = -dtp. (so) 

r 

The volume element is given by 

dV = rdrdOdz = {r/R)drdsdz, (81) 


where ds = RdO. We then integrate both hand sides of 
(l80t by {r/R)dr in a region R — ej^ ^ R + ej^ and then 
we obtain two dimensional conservation law. 

dsJlot + = -5tP, J P dsdz = 1, (82) 


where 


JL ^ 


f'R-\-e/2 


rJ^ dr, 


lR-e/2 
pR+e-12 

Jl, = / rr/R dr, 

jR-e/2 
rR+e/2 

P= rp/R dr, 

JR-e/2 


(83) 

(84) 

(85) 


where the boundary condition J'' = 0atr = i?±e/2 and 
ds = RdO are utilized. We have 


7® — 


2mi 


f'R-V^ j 2 


'R-t/2 


dr 


{'ip*dei> - i’dsip*)- 


( 86 ) 


FIG. 2: Bending ribbon (side view) with thickness e. 

The physical space is inner ribbon with 0 < z < L, and 
—e/2 < r — R < -|-e/2. Our starting equation is an usual 


three dimensional Schrddinger equation (I20p. 
obtain the conservation law. 

Then we 

“ dt^ 

(73) 


(74) 


The separation of variable (l24l) , ([26|) gives 

ip{r,e,z) = {g/Gy/'^p{e,z)xir) = \f^‘p{e,z)x{r), 

V r 

(87) 

where 

J dr Ixp = 1, J RdOdz \p\^ = J dsdz \ip\^ = 1, 

ds( 2 ) = R^dO^ + dz"^, gij = diag(i?^, 1), y/g = R. 
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Then we obtain 

= -,\xir)?i^*ds^-^ds^l. ( 88 ) 

Equation (l54)) gives the explicit form of x- By using 
C = 9 °/e = (r - R)/e = i?a;/e, 

we obtain 


This is equal to in (1711) , and we have no geometric flow 
in this straight direction. Geometric flow is contained as 
the part of the integrand in equation (IMll : the second 
order of expansion We can transcribe (1551) into the 
form 

JL = ^ < ( 7 )' > (93) 

f‘R-\-ej2 

<f>qm= / /(r) |x(r)pdr. (94) 

jR-el2 


— 


HR 


('-\-ej2R 


dx 


cos^(( 7 V 7 ri?/e)x), 


mieJ_^/2R (l + a;)2 

X {(p*ds(p- (pdsif*), 


(89) 


for odd N. For even N we just change cos to sin. 
We expand the integrand as 


5. Conclusion 

We have discussed the conservation law in an effective 
two-dimensional system between two curved surfaces T,' 
and E separated by a small distance e. We found that the 
anomalous flow depends on the curvature of the surface 
E. In the classical diffusion process, we have 


1 


(l+x)2 

Then we obtain 


— 1 — 2 x “h 3x 


(90) 


JL = + (91) 

Imi 


•Jn — ~ 


e^D 

IT 


(3k*"^k7 - 2kC^) 


dqi 


1 

-Q ^ - 

2^ dq^ 


as shown in [^. 

In the quantum process we instead obtain 


(95) 


for both of even and odd N. If we use coordinate g* = 
( 5 , 2 ), we have Kss = 1/7?, Hsz = Kzz = 0. 

Then we find this is completely equal to J® -I- Jq appeared 
in m and dZH). 

On the other hand, 

JLt = dr \x{r)\'^){(p*d^(p-(pd^(f*) 

2mi 

= (92) 


dh = ^(1 - ^ 7 ^) - 7^9jT)- (96) 

The classical anomalous flow and quantum mechanical 
anomalous flow are somewhat similar. Both start from 
and are proportional to the gradient of the field 
except the last term in the classical flow. The recent 
nano-technology made it possible to fabricate complicate 
devices. Then this kind of anomalous flow might play an 
important role in such a “geometrical” device. 


[1] J. Faraudo, J. Chem. Phys, 116 (2002) 5831-5841; J. 
Balakrishnan, arXiv:physics/0308089, 25 Aug 2003; N. S. 
Gov, Phys. Rev. E 73 (2006) 041918; A. Naji, F. Brown, 
J. Chem. Phys. 126 (2007) 235103; E. Reister, U. Seifert, 
arXive:cond-mat/0503568, 23 Mar 2005; C. Venkatara- 
man, T. Sekimura, E. Gaffney, P. Maini, A. Madzvamuse, 
Phys. Rev. E 84 (2011) 041923; R. C. Priego, P. C. Villar¬ 
real, S. E. Jimenez, J. M. Alcaraz, arXivexond-mat.stat- 
mech/1211.5799v2 11 Jun 2013; P. C. Villarreal, A. V. 
Balbuena, J. M. Alcaraz, R. C. Priego, S. E. Alvarez, J. 
Chem. Phys. 140 (2014) 214115; T. Balois, C. Chatelain, 
M. B. Amar, J. R. Soc. Interface 11 (2014) 97; T. Do- 
browolski, Eur. Phys. J. B 86 (2013) 346; R. G. Morris, 
Phys. Rev. E 89 (2013) 062704. 

[2] N. Ogawa, Phys. Rev. E81 (2010) 61113. 

[3] N. Ogawa, Phys. Lett. A377 (2013) 2465-2471. 

[4] R. C. T. da Costa, Phys. Rev. 23 (1981) 1982; J. Tolar, 


1988 Lecture Notes in Physics 313, ed. H. D. Doever, J. D. 
Henning, T. D. Raev (Springer-Verlag, Berlin, Heidelberg) 
268. 

[5] N. Ogawa, K. Fujii, K. P. Kobushkin, Prog. Theor. Phys. 
83 (1990) 894; N. Ogawa, K. Fujii, N. M. Chepilko, K. P. 
Kobushkin, Prog. Theor. Phys. 85 (1991) 1189; N. Ogawa, 
Prog. Theor. Phys. 87 (1992) 513. 

[6] K. Fujii, N. Ogawa, Prog. Theor. Phys. 89 (1993) 575. 

[7] Gorria C, Gaididei YuB, Soerensen MP, Christiansen PL, 
Caputo JG. Phys. Rev. B69 (2004) 134506; T. Dobrowol- 
ski. Archives of Materials Science and Engineering 39 2 
(2009) pll6-120; T. Dobrowolski, Phys. Rev. E79 (2009) 
046601; T. Dobrowolski, Canadian J. Phys. NRG, 88 9 
(2010) p627-633; T. Dobrowolski, Discrete and Continu¬ 
ous Dynamical Systems - Series S 4 5 (2011) pl095-1105; 
T. Dobrowolski, Eur. Phys. J. B 86 (2013) 346. 








